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ABSTRACT 

Recent observations suggest that gamma ray bursts (GRBs) and their 
afterglows are produced in star formation regions in distant galaxies by highly 
relativistic jets that happen to point in our direction. Relativistic beaming 
collimates the emission from the highly relativistic jets into small solid angles 
along the jet direction. It implies that we are seeing only a small fraction of 
the events that produce GRBs. The observed GRB rate then requires an event 
rate which is comparable to the birth rate of neutron stars (NS). The highly 
relativistic jets sweep up ambient matter along their trajectories, accelerate it 
to cosmic ray (CR) energies and disperse it in hot spots which they form when 
they stop in the galactic halo. With an event rate comparable to the NS birth 
rate, such events in our Galaxy may be the main source of Galactic cosmic 
rays at all energies. Internal interactions and/or external interactions of these 
jets with high column density matter and/or radiation at their production sites 
or along their trajectories can produce high energy gamma rays and neutrinos 
that are highly beamed along the jet direction. Jetted GRBs, like blazars, 
may be much more fluent in high energy gamma rays and neutrinos than in 
MeV gamma rays. But, TeV gamma rays from large cosmological distances are 
unobservable because of their attenuation by electron-positron pair creation on 
the intergalactic infrared background radiation. However, high energy neutrinos 
from distant GRBs may be observed with large surf ace/ volume telescopes which 
are under construction. TeV gamma rays and high energy neutrinos may also be 
detected from relatively nearby GRBs by the existing moderate size detectors, 
but with a much smaller rate. 
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1. The Energy Crisis Of Spherical GRBs 

Thanks to the precise and prompt localization by the Italian-Dutch satellite, BeppoSAX 
(see, e.g., Costa et al. 1997), long lived GRB afterglows spanning the wavelength range 
from X-ray to radio have now been detected in more than a dozen GRBs. They led to 
the redshift measurements, z=0.69, 0.835, 3.42, 1.096, 0.966, 1.61, 1.62 of GRBs 970228 
(Kulkarni et al. 1999) 970508 (Metzger et al. 1997), 971214 (Kulkarni et al 1998), 980329, 
980613 (Djorgovski et al 1999), 980703 (Djorgovski et al. 1998), 990123 (Anderson et al. 
1999; Kulkarni et al. 1999), 990510 (Vreeswijk et al. 1999), respectively, from absorption 
lines in their optical afterglows and/or emission lines from their host galaxies. In addition, 
strong suppression has been observed with the Hubble Space Telescope in the spectrum 
of the host galaxies of GRBs 970228 and 980329 at wavelengths below 700 nm. If it is 
due to absorption in the Lya forest (Fruchter 1999), then their redshifts are near z~ 5. 
These measured/estimated redshifts indicate that most GRBs take place at very large 
cosmological distances. For instance, assuming a zero cosmological constant (f^ = 0), the 
luminosity distance 

c 2[2 - Q M (1 - z) - (2 - n M )Vl + M z] m 

D L=u 7T2 > U) 



h n 2 M 

with the present canonical values for the cosmological parameters, H = 65 km cm -1 Mpc~ 4 
and f2 M = 0.2 (which will be assumed in this paper), yield D L ~ 5 x 10 28 cm for z=2. 
The typical observed GRB ffuence (F 7 ~ 10" 5 erg cm -2 ) and their large distances imply 
enormous energy release in gamma rays, 

4vrD 2 F 7 „ 
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if their energy release is isotropic as used to be assumed/advocated by the standard 
fireball models of GRBs and GRB afterglows (e.g., Piran 1999 and references therein). In 
particular, the large fluence F 7 5.1 x 10~ 4 erg cm -2 (Kippen et al. 1999) and redshift 
z=1.61 of GRB 990123 yield E 7 ps 3.4 x 10 54 erg. Such enormous energy release in gamma 
rays alone, implies an "energy crisis" for spherical GRBs (Dar 1998): The short duration 
and the very large energy release in GRBs indicate that they are powered by gravitational 
collapse of compact stars. But, the energy release in such events falls short of that required 
to power GRBs like 971214, 980329 and 990123, 990510 if they were isotropic. Furthermore, 
all the known luminous sources of gamma rays (quasars, radio galaxies, active galactic 
nuclei, accreting binaries, pulsars, supernova explosions, supernova remnants) exhibit rather 
a modest efficiency, rj < 10~ 4 , in converting gravitational, kinetic or thermonuclear energy 
into gamma rays. If GRBs have a similar efficiency for converting the energy release from 
their central engine to gamma rays, then the energy crisis is common to most GRBs. 
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2. Jetted GRBs 

2.1. No Energy Crisis For Jetted GRBs 

Various authors have pointed out that the energy crisis in isotropically emitting GRBs 
is avoided if GRBs are beamed into a small solid angle, Ail <C Att such that their total 
energy release in gamma rays is 

Beaming of gamma rays from GRBs is possible if the highly relativistic ejecta (Lorentz 
factor T = 1 / y/1 — (3 2 ^> 1) that produces the GRB is beamed into a cone (conical beaming) 
of solid angle Ail <C 4ir, or if the ejecta is jetted - namely, if after initial expansion 
the ejected cloud/plasmoid maintains nearly a constant cross section. Conical beaming 
(e.g., Mochkovich 1993; Rhoads 1997) can solve the " energy crisis" of GRBs by reducing 
their inferred energies by the ratio Af2/47r <C 1. But, it suffers from other deficiencies of 
isotropically emitting GRBs (Dar 1998). Jetting the ejecta (e.g., Shaviv and Dar 1995, 
Dar et al. 1998, Dar 1998) solves the energy crisis and can also explain the short time 
variability of GRB light curves, the versatility of their afterglows, the absence of a simple 
scaling between them and their sudden decline in some GRBs (970508, 990123, 990510): 

The emission from a highly relativistic plasmoid which is isotropic in its rest frame 
and has a power-law spectral shape, F u i = Au'~ a , is collimated in the lab frame to small 
emission angles 9 ~ 1/T relative to its direction of motion, according to 

Fu = (1 + r 2^ 2 )3 F ^'=Ki+r 2 e 2 )/2r • (4) 

Thus, the observed flux from a plasmoid with a typical spectral index a ~ 0.7 that moves 
with a Lorentz factor T ~ 10 3 at an angle 9 < 1/T relative to the line of sight is amplified 
by approximately T 3+a ~ 10 11 . This amplification within a solid angle Ail ~ ir/T 2 can 
explain why highly relativistic jets with bulk motion Lorentz factors T ~ 10 3 , total kinetic 
energy Ek ~ 10 52 erg, and conversion efficiency i] > 10~ 4 into gamma rays can produce 
GRBs with equivalent isotropic energy of E 7 = iiE^Tr/Ail > 4 x 10 54 erg, as observed for 
GRB 990123 (Kippen et al. 1999). 



2.2. The Beaming Angle Of GRBs 



The enormous release of energy in GRBs during a short time suggests that they are 
energized by collapse of compact stars (Blinnikov 1984, Paczynski 1986) due to mass 
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accretion (Goodman, Dar and Nussinov 1987; Dar et al. 1992) or phase transition (e.g., 
Dar 1999a). If GRBs are produced, e.g., by gravitational collapse of neutron stars (NS) 
to quark stars (QS) when they cooled and spun down sufficiently (e.g., Dar 1999a; Dar 
and De Riijula 1999), then the GRB rate is comparable to the NS birth-rate. The NS 
birth-rate is estimated to be Rns ~ 0.02 y" 1 in Milky Way like galaxies (van den Bergh 
and Tamman 1991). From the observed rate of GRBs, Rqrb[UNIV] ~ 10 3 y _1 in the whole 
Universe, it was estimated that the rate of observable GRBs in Milky Way like galaxies is 
R G rb[MW] ~ 1(T 8 y- 1 (e.g., Wijers et al. 1997). The beaming angle of GRBs therefore 
must satisfy 

R G rb[MW] « 2(Afi/47r)R NS [MW] (5) 

where we assumed that two opposite jets are ejected in every NS collapse. Hence, 
Af2 ~ tt/T 2 ~ 7r x 1CT 6 . It implies that their bulk motion Lorentz factor is T ~ 10 3 . Such 
values have been inferred also from the absence of a break due to 77 — > e + e~ in GRB 
spectra (e.g., Baring and Harding 1997), from the peak energy of GRBs and from GRB 
duration and substructure (e.g., Shaviv and Dar 1995). Such strong beaming implies that 
we observe only a very small fraction, ~ 10~ 6 , of the events that produce GRBs. I will call 
these events cosmological GRBs (CGRBs) if they occur in distant galaxies and "Galactic" 
GRBs (GGRBs) if they occur in our Milky Way (MW) galaxy. 



2.3. The Jet Energy 

Consider gravitational collapse that leads to the birth of a pulsar (e.g., gravitational 
collapse of NS to QS due to a phase transition of cold and highly compressed neutron 
matter to Bose condensate of diquark pairs [Dar 1999a, Dar and De Riijula 1999] or the 
birth of a pulsar in a supernova explosion [Cen 1999]) and perhaps to the ejection of two 
opposite highly relativistic jets. If momentum imbalance in the ejection of the relativistic 
jets (and not asymmetric neutrino emission) is responsible for the observed large mean 
velocity (Lyne and Lorimer 1994), v as 450 ± 90 km s -1 , of slowly spinning pulsars, then 
momentum conservation implies that the difference in the kinetic energy of the jets satisfies 

AE jct > cP ns ~ vM NS c ~4x 10 51 erg, (6) 

where we used the typical observed mass of NSs, M^s ~ 1.4M Q . If AEj Ct <C Ej Ct , then the 
kinetic energy of the jets must be Ej et ~ 10 52 erg or larger. If T rs 10 3 then the ejected jet 
(plasmoid) has a mass M jct ~ 1.5 x 10~ 6 M NS ~ 2.1 x 1O~ 6 M ~ 0.7M Ea rth- Even if only a 
fraction rj ~ 10~ 4 of the jet kinetic energy is radiated in 7-rays, the inferred "isotropic" 
7-ray emission in GRBs is E isot ~ Ar]T 2 Ej Ct ~ 4 x 10 54 erg, while the true 7-ray emission is 
only E 7 ~ 10 48 erg. 
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2.4. Jet Formation 

Relativistic jets seem to be emitted by all astrophysical systems where mass is 
accreted at a high rate from a disk onto a central compact object. Highly relativistic jets 
were observed in galactic superluminal sources, such as the microquasars GRS 1915+105 
(Mirabel and Rodriguez 1994; Mirabel and Rodriguez 1999) and GRO J 165-40 (Tingay et 
al. 1995) where mass is accreted onto a stellar black hole (BH), and in many active galactic 
nuclei (AGN), where mass is accreted onto a supermassive BH. Mildly relativistic jets from 
mass accretion are seen both in AGN and in star binaries containing NS such as SS433 
(e.g., Hjellming and Johnston 1988). The ejection of highly relativistic jets from accreting 
or collapsing compact stellar objects is not well understood. Therefore their properties 
must be inferred directly from observations and/or general considerations: High-resolution 
radio observation resolved the narrowly collimated relativistic jets of microquasars into 
blobs of plasma (plasmoids) that are emitted in injection episodes which are correlated with 
sudden removal of the accretion disk material (Rodriguez and Mirabel 1998). After initial 
expansion, these plasmoids seem to retain a constant radius of R p ~ 10~ 3 pc. The emission of 
Doppler shifted Hydrogen Lya and Iron Ka lines from the relativistic jets of SS433 suggest 
that the jets are made predominantly of normal hadronic plasma. Moreover, simultaneous 
VLA radio observations and X-ray observations of the microquasar GRS 1915+105 indicate 
that the jet ejection episodes are correlated with sudden removal of accretion disk material 
into the relativistic jets (Mirabel and Rodriguez 1999). Highly relativistic jets, probably, 
are also ejected in the birth or collapse of NSs due to mass accretion or phase transition. 
But, because the accretion rates and magnetic fields involved are much larger compared 
with those in quasars and microquasars, the bulk motion Lorentz factors of these jets 
may be much higher, perhaps, T ~ 10 3 as implied by the above consideration and GRB 
observations. When these highly relativistic jets happen to point in our direction they 
produce the observed cosmological GRBs and their afterglows. They look like the Galactic 
micro copies of blazer ejections and therefore will be called "microblazers" . In fact, when 
the light curves and energy spectra of blazar flares and of microquasar plasmoids are scaled 
according to the Lorentz factors expected for GRB ejecta, they look quite similar to GRBs 
(Dar 1999b). 

The high collimation of relativistic jets over huge distances (up to tens of pc in 
microquasars and up to hundreds of kpc in AGN), the confinement of their highly 
relativistic particles, their emitted radiations and observed polarizations, all indicate that 
the jets are highly magnetized, probably with a strong helical magnetic field along their 
axis. Magnetic fields as strong as a few tens mili Gauss in the jet rest frame have been 
inferred from microquasar observations (Mirabel and Rodriguez 1999), while hundreds 
Gauss were inferred for GRB ejecta (assuming equipartition of energy between internal 
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kinetic and magnetic energy). The UV light and the X-rays from the jets ionize the ISM in 
front of them. The jet material and the swept-up ionized ISM material in front of the jet 
can be accelerated by the Fermi mechanism to a power-law energy distribution that extends 
to very high energies, as inferred from the observed radiations from jets. The interactions 
of these high energy particles in the jet and/or their interactions with the external medium, 
produce the GRBs and their afterglows: 



2.5. Jet Production of Gamma Rays 

The GRB jets may consist of pure e + e~ plasmoids or of normal hadronic gas/plasma 
clouds. The GRB can be produced by electron synchrotron emission. If the jet consist of 
a single plasmoid, then individual 7-ray pulses that combine to form the GRB light curve 
can be produced by internal instabilities or by interaction with inhomogeneous external 
medium. If the jets consist of multiple ejections of plasmoids, then the GRB pulses may 
be produced when later ejected plasmoids collide with earlier ejected plasmoids that have 
slowed down by sweeping up the interstellar medium in front of them. But, such scenarios 
do not seem to provide a simple explanation why the GRB emission is peaked near ~ 
MeV photon energy. Other GRB emission mechanisms, however, can provide such an 
explanation: 

If the highly relativistic plasmoid consists of a pure e + e~ plasma, then inverse Compton 
scattering of stellar light (hu = e cv x 1 eV) by the plasmoid can explain the observed 
typical 7 energy (e 7 ~ 4r 2 e eV /3(l + z) MeV), GRB duration (T ~ R SFR /2cr 2 ~ 50 s), 
pulse duration (t p ~ R p /2cr 2 ~ 150 ms), fluence (F 7 ~ 10~ 5 erg cm -2 ), light curve and 
spectral evolution of GRBs (Shaviv and Dar 1995; Shaviv 1996; Dar 1998). For instance, 

_ q T Ne 7 E jet (l + z) _ lQ- 5 z 2 N 22 r 3 e cv E 52 erg 
7 ~ rm c c 2 D 2 Afi ~ D 2 9 cm 2 1 ' 

where D = D 29 x 10 29 cm is the luminosity distance of the GRB at redshift z, z 2 = (1 + z)/2, 
N = N 22 x 10 22 cm" 2 is the column density of photons along the jet trajectory in the star 
forming region, a T = 0.65 x 10~ 24 cm -2 is the Thomson cross section, Ej et = E 52 x 10 52 erg 
and T = T 3 x 10 3 . 

If the plasmoid consists of normal crust material of neutron stars (Doppler-shifted 
K a iron line was detected from the jets of SS433), then photoabsorption of stellar 
light by partially ionized heavy metals like iron (Doppler-shifted to X-rays in the jet 
rest frame) and its reemission as 7 rays (iron X-rays lines in the jet rest frame) yield 
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e 7 ~ Te x /(1 + z) ~ MeV) in the observer frame and 

^ g a Ne 7 E jet (l + z) ^ 10~ 5 Z2ai9N22exr 3 E 5 2 erg 
7 ~ rM Fe c 2 D 2 Afi ~ D 2 9 cm 2 1 J 

where <r a = cr 19 x 10~ 19 cm 2 is the mean photoabsorption cross section of X-rays by partially 
ionized iron. 



3. GRB Afterglows 

The afterglows of GRBs may be synchrotron emission from the decelerating plasmoids 
(e.g., Chiang and Dermer 1997), and then they are highly beamed and may exhibit 
superluminal velocities (c < vj_ < Tc) during and right after the GRB. The deceleration of 
a mildly relativistic spherical ejecta from NS collapse to QS may also produce spherical 
supernova like light curve (many planetary nebulae, e.g., NGC 7009, NGC 6826, and some 
SNRs (including, perhaps, SNR 1987A) show antiparallel jets superimposed on a spherical 
explosion). In the rest frame of the decelerating plasmoid, the synchrotron spectra can be 
modeled by convolving the typical electron energy spectrum (E~ p at low energies up to 
some "break energy" where it steepens to E~ p_1 and cuts off exponentially at some higher 
energy due to synchrotron losses in magnetic acceleration) with the synchrotron Green's 
function (see, e.g., Meisenheimer et al. 1989). In the observer frame it yields spectral 
intensity (Dar 1998) 

T ,,a + /3 , -0.75±0.25 + -1.25±0.08 / n \ 

where a = — (p — l)/2 and f3 = — (p + 5)/6 and where I assumed p = 2.5 ± 0.5 for Magnetic 
Fermi acceleration. This prediction is in agreement with observations of GRB afterglows. 
Moreover, the glows of microquasar plasmoids and radio quasar jets after ejection, and of 
blazar jets after flares, show the same universal behavior as observed in GRBs afterglows. 
For instance, the glows of the ejected plasmoids from GRS 1915+105 on April 16, 1994 near 
the source had a = —0.8 ± 0.1 and (3 = —1.3 ± 0.2 (Rodriguez and Mirabel 1998), identical 
to those observed for SS 433 (Hjellming & Johnston 1988) and for the inner regions of jets of 
some radio galaxies (e.g., Bridle & Perley 1984). When the jets spread, the spectral index of 
their power-law time decline changed to f3' = 2.6 ± 0.4 (e.g., Mirabel and Rodriguez 1999). 
Thus, their overall time decline can be described approximately by I ~ t _/3 /[l + (t/to) /3 '~ /3 ] 
where to is the time when the jet begins to spread. Indeed, such a behavior has been 
observed also in the afterglows of some GRBs (e.g., GRB 990510; Vreeswijk et al 1999). 
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4. Galactic GRBs - The Main Source of Cosmic Rays? 

According to the current paradigm of cosmic rays (CRs) origin, CR nuclei with energies 
below 3 ■ 10 15 eV (the "knee") are accelerated in Galactic SNRs (Ginzburg 1957) and those 
above 3 x 10 18 eV (the "ankle"), for which a disk origin is unlikely due to their isotropy, in 
sources far beyond our Galaxy (Burbidge 1962). However, recent observations suggest that, 
perhaps, SNRs are not the main source of Galactic cosmic rays and that the CRs above the 
ankle are not extragalactic: 

- Measurable fluxes of high energy gamma rays from interactions of cosmic ray nuclei 
in SNRs, as expected in models of SNR acceleration of CRs, were not detected from nearby 
SNRs (Prosch et al. 1996; Hess et al. 1997; Buckley et al. 1998). 

- The expected galactocentric gradient in the distribution of high energy gamma rays 
(> 100 MeV) from interactions of CRs from SNRs in the Galactic interstellar medium is 
significantly larger than observed by the EGRET detector on board CGRO (Hunter et al. 
1997; Strong and Moskalenko 1998). 

- Diffusive propagation of CRs from the observed Galactic distribution of SNRs yield 
anisotropies in the distribution of CRs above 100 TeV in excess of the observed value 
(Aglietta et al. 1995) by more than an order of magnitude (Ptuskin et al. 1997). 

- The absence (Takeda 1998) of the "GZK cutoff" in the intensity of CRs at energies 
above ~ 10 20 eV due to interactions with the cosmic microwave radiation (Greisen 1996; 
Zatsepin & Kuz'min 1966) have brought into question (e.g., Hillas 1998) their hypothesized 
extragalactic origin. 

Relativistic jets are efficient CRs accelerators (e.g., Mannheim and Biermann 1992, Dar 
1998b). A modest fraction of the total energy injected into the MW in jets from Galactic 
GRBs at a rate similar to the NS birth/collapse rate in the MW, that is converted to CR 
energy, can supply, ~ 1.5 x 10 41 erg s -1 , the estimated Galactic luminosity in CRs (Drury 
et al., 1989). Thus, Dar and Plaga (1999) have recently proposed that Galactic GRBs 
(GGRBs) are the main source of the CRs at all energies and consequentlty no GZK cutoff 
is expected in the CR spectrum: 

The highly relativistic narrowly collimated jets/plasmoids from the birth or collapse of 
NSs in the disk of our Galaxy that are emitted with Ej Ct ~ 10 52 erg perpendicular to the 
Galactic disk, stop only in the Galactic halo, when the rest mass energy of the swept-up 
ambient material becomes > their initial kinetic energy. Through the Fermi mechanism, 
they accelerate the swept-up ambient matter to CR energies and disperse it into the halo 
from the hot spots which they form when they finally stop in the Galactic halo (Fig.l). 



- 9- 




Fig. 1. — A highly schematic sketch of the Dar-Plaga paradigm. The birth or collapse of NSs 
in the disk of our Galaxy leads to an ejection of two opposite jets that produce "hot spots" 
when they stop in an extended Galactic halo. 

The typical equipartition magnetic fields in such hot spots may reach 
B ~ (SEjct/Rp) 1 / 2 ~ 1 G. Synchrotron losses cut off Fermi acceleration of CR nuclei 
with mass number A at E~ rA 2 Z~ 3 / 2 (B/G)- 1/2 x 10 20 eV. Particle-escape cuts 
off Fermi acceleration when the Larmor radius of the accelerated particles in the 
plasmoid rest frame becomes comparable to the radius of the plasmoid, i.e., above 
E ~ rZ(B/G)(Rp/0.1 pc) x 10 20 eV. Consequently CR with E > Z x 10 20 eV can no longer 
be isotropized by acceleration or deflection in hot spots with T ~ 1. 

Fermi acceleration in the highly relativistic jets from GRBs (T ~ 10 3 )can produce 
a broken power-law spectrum, dn/dE ~ E~ a , with f3 ~ 2.2 below a knee around 
Eknee ~ A PeV and (3 ~ 2.5 above this energy (Dar 1998b). Spectral indices a ~ 2.2 were 
obtained also in numerical simulation of relativistic shock acceleration (e.g., Bednarz and 
Ostrowski 1998) . Galactic magnetic confinement increases the density of Galactic CRs 
by the ratio CTh/Rc where 7h(E) is the mean residence time in the halo of Galactic CRs 
with energy E, and Rg ~ 50 kpc is the radius of the Galactic magnetic-confinement region. 
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With the standard choice for the energy dependence of the diffusion constant (observed, 
e.g., in solar-system plasmas) one gets: r h oc (E/Z)~°- 5 . Consequently, the energy spectrum 
of CRs is predicted to be 

dn/dE ~ C(E/E kncc )- p (fO) 

with p ~ a + 0.5 ~ 2.7 (~ 3) below (above) the knee. This power-law continues as long as 
the Galactic magnetic field confines the CRs. 

Part of the kinetic energy released by GGRBs is transported into the Galactic 
halo by the jets. Assuming equipartition of this energy, without large losses, between 
CR, gas and magnetic fields in the halo during the residence time of CR there, the 
magnetic field strength Bh in the halo is expected to be comparable to that of the disk 
Bh ~ (2Lmw[CR]ti 1 /R^) 1 / 2 ~ 3 /xG where ~ 5 x 10 9 y is the mean residence time of the 
bulk of the CRs in the Galactic halo. Cosmic rays with Larmor radius larger than the 
coherence length A of the halo magnetic fields, i.e., with energy above 

E ank ie ~ 3 x 10 18 (ZB h /3/iG)(A/kpc) eV, (11) 

escape Galactic trapping. Thus, the CR ankle is explained as the energy where the mean 
residence time Th(E) of CRs becomes comparable to the free escape time from the halo 
Tf r ee ~ 1.6(i?h/50 kpc) x 10 5 years. Therefore, the spectrum of CRs with energies above the 
ankle, that do not suffer Galactic magnetic trapping, is the CRs spectrum produced by the 
jet, i.e., 

dn/dE ~ C(E ank i e /E knee )~ 3 (E/E ank i c )~ 2 ' 5 ; E > E ank i c . (12) 
Eqs. 10-12 describe well the overall CR energy spectrum. 



5. High Energy Gamma Rays 

The observed radiations from blazars, microquasars and GRBs indicate that their 
highly relativistic jets contain high energy charged particles with a power-law energy 
distribution dn p /dE ~ AE~ Q with a ~ 2.2, that extends to very high energies. Such 
distributions can be formed in the highly relativistic jets through Fermi acceleration of 
swept up material. This power-law distribution of energetic protons is boosted and beamed 
into a solid angle AQ tv tt /T 2 in the lab frame. GRB afterglows suggest that GRBs are 
produced in star formation regions, probably molecular clouds. The typical column density 
of such clouds is N p = 10 24±1 cm -2 . The clouds must be highly ionized along the line of sight 
to the GRB by the enormous fluxes of beamed X-rays and UV radiations from the GRB. 
Interaction of the highly relativistic GRB jets with the high column density of this ionized 
gas (or with diffuse matter at their production sites) can produce high energy gamma rays 
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Fig. 2. — The observed flux of cosmic rays (thin line) as a function of primary energy E is 
well described by a power law that changes its slope sharply at only two energies the "knee" 
and the "ankle" . At energies below the ankle it is enhanced (by a factor (E/E ank i e )~ - 5 ) over 
the GGRB source spectrum (thick line, a power law with differential power law index of -2.2 
below the knee and ~ -2.5 above it) by way of trapping in the Galactic halo magnetic fields. 



through pp — > 7r°(?7 )X; n Q (rf) — > 27. The cross section for inclusive production of high 
energy 7-rays with a small transverse momentum, cpx = Et < 1 GeV in pp collisions (e.g., 
Neuhoffer et al. 1971; Boggild and Ferbel 1974; Ferbel and Molzon 1984) is well represented 
by 

— H^ET « (l/27rp T )e- ET / E ° f(x), (13) 
<7 in d 2 p T dE 7 

where E is the incident proton energy, <r in ps 35 mb is the pp total inelastic cross section at 
TeV energies, E ~ 0.16 GeV and f(x) ~ (1 — x) 3 /^^ is a function only of the Feynman 
variable x = E 7 /E, and not of the separate values of the energies of the incident proton 
and the produced 7-ray (Feynman scaling). The exponential dependence on E T beams 
the 7-ray production into < Et/E ~ 0.17/T along the incident proton direction. When 
integrated over transverse momentum the inclusive cross section becomes <T i ~ 1 dcr/dx ps f(x). 
If the incident protons have a power-law energy spectrum, dn p /dE ps AE~ a , then, because 
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of Feynman scaling, the produced 7-rays have the same power-law spectrum: 

dm 



dE 7 ~ NpC7in 7e 7 dE dE 7 



<lll|> ' (IK^.X^^p-AK,-. 



(14) 



where N p is the column density of the target and g P7 = Jq 1 x a_1 f(x)dx ks 0.092 for a ~ 2.2. 
Consequently, the collimated flux of high energy gamma rays produced by a GRB jet with 
initial kinetic energy E = E52 x 10 52 erg that propagates through a molecular cloud of 
typical column density N p = N23 x 10 23 cm" 2 is given by 



dm 6 x 10~ 6 E 

r^i 

dE 



5 2 N 2 3r| 



n 2 



;i + z) 



2-Q 



E 



TeV 



-2.2 



e -r(z,E) cm -2 Te y-1 



(15) 



where D(z) = D 29 x 10 29 cm, is the luminosity distance to the GRB and r(z, E) is the 
optical depth to the GRB (redshift z) at energy E. The fluxes of high energy gamma rays 
from GRBs at z ~ 2 are attenuated strongly (r > 1) for E> 20 GEV. For not very distant 
GRBs, e.g., z< 0.5, the gamma ray flux is not attenuated strongly at energy E < 100 GeV 
(Salomon and Stecker, 1998). GRBs with z< 0.1 (D L < 500Mpc) can be visible in TeV 
gamma rays. But, their expected rate is only 



R G rb(z < 0.1) ~ R GRB [L*KV c (z < 0.1) ~ 0.1 y 



-1 



(16) 



where Rqrb[L*] ~ 2 x 10~ 8 y _1 is the estimated mean rate of GRBs in L* 1O 1O L galaxy 
for z< 0.1, p L ~ 1.2 x 1O 8 L Mpc~ 3 is the luminosity density in the local Universe and 
V c (z < 0.1) ~ 5 x 10 8 Mpc 3 is the comoving volume within z< 0.1. GeV gamma rays from 3 
very luminous GRBs have been reported by the EGRET detector on board CGRO (Dingus 
et al 1994; Dingus 1995), consistent with the above predictions. However, they could have 
also been produced by inverse Compton scattering of GRB photons from energetic electrons 
in the GRB jets. Only the detection of high energy neutrinos from GRBs can establish the 
hadronic production origin of high energy photons from GRBs, i.e., the hadronic nature of 
the GRB jets. 



6. High Energy Neutrinos From GRBs 

Hadronic production of photons in diffuse targets is accompanied by neutrino 
emission mainly through hadronic production of mesons that decay into neutrinos, e.g., 
pp — > — > (J^Vfj, ; fi^ — > e ± h , f j i u e . Analytical calculations (e.g., Dar 1983, Dar 1984, 
Lipari 1993) show that a proton power-law spectrum, dn p /dE = AE~ a with a power 
index a ~ 2.2, generates power-law spectra of 7-rays and v^s that satisfy approximately, 
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dn„/dE « 0.80dn 7 /dE (Dar and Shaviv, 1996). Consequently, 



dn, 5 x 10- 6 E, 2 No,r 



2 



52^23*- 3/-, , „\2-a 



dE Dl 9 

Thus, we predict that the high energy 7-ray emission from GRBs is accompanied by 
emission of high energy neutrinos with similar fluxes, light curves and energy spectra. The 
expected number of events from a GRB in a deep underwater /ice telescope is 

Nevents « SN A J R^^dE.dE (18) 

where S is the effective surface area of the telescope, Na is Avogadro's number, a Vil is the 
inclusive cross section for z/^p — > /xX, and R M is the range (in gm cm -2 ) of muons with 
energy in water/ice. The number of events is not sensitive to the detector energy 
threshold, if it is below 300 GeV where both the muon range and production cross section 
increase linearly with energy and yield a detection probability in ice/water, which increases 
like ~ 10 _6 (E/TeV) 2 below 300 GeV. Using the neutrino cross sections that were calculated 
by Gandhi et al (1998) and neglecting detector threshold effects and neutrino attenuation 
in Earth (which becomes important only above 100 TeV), we predict that the number of 
neutrino events in deep underground ice/water neutrino telescopes per 1 km 2 is 

N events « 1.3 x (1 + z) -o.2^3T| km _ 2 (ig) 

Thus, a relatively nearby GRB (z=0.5) may generate ~ 40E 52 N 2 3 upgoing muon 
events in underwater /ice telescope per 1 km 2 area and only ~ IE52N23 events if it is 
at z=2. The expected time length of these neutrino bursts from GRBs is typically, 
t R c i/cr 2 ~ Rio x 10 3 s, where R c i = IOR10 P c is the size of the molecular cloud. Such 
events can be distinguished from the atmospheric neutrino background by their directional 
and time coincidence with the GRBs and establish the hadronic nature of the relativistic 
jets from GRBs. 

Unlike the neutrino bursts from nearby supernova explosions, the arrival times of 
v's from GRBs, which are spread over t > t3 x 10 3 s, yield only poor limits on neutrino 



1 + z)' 



E 



TeV 



-2.2 



cm" 2 TeV" 1 . (17) 



masses and lifetimes: im,c 2 > ^/2t/TE„ ~ y t 3 /Ti [E„/Tev] x 10 5 eV, where Ti is the 
GRB lookback time in units of lOGy. This limit cannot compete with the cosmological 
limit, I] m n U c 2 < 94f2 M h 2 eV 8 eV, for long lived neutrinos. The neutrino arrival times 
can be used, however, to improve the limit from Supernova 1987A (LoSecco 1987) on the 
equivalence principle of General Relativity. 



Results presented here are based on an ongoing collaboration with A. De Rujula and 
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R. Plaga. This research was supported by the fund for the promotion of research at the 
Technion. 
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